home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SGESL.z / SGESL
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSGGGGEEEESSSSLLLL((((3333FFFF))))                                                            SSSSGGGGEEEESSSSLLLL((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SGESL   - SGESL solves the real system A * X = B  or  TRANS(A) * X = B
  10.      using the factors computed by SGECO or SGEFA.
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE SGESL(A,LDA,N,IPVT,B,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAA REAL(LDA, N)
  20.         the output from SGECO or SGEFA.
  21.  
  22.      LLLLDDDDAAAA INTEGER
  23.         the leading dimension of the array  A .
  24.  
  25.      NNNN INTEGER
  26.         the order of the matrix  A .
  27.  
  28.      IIIIPPPPVVVVTTTT INTEGER(N)
  29.         the pivot vector from SGECO or SGEFA.
  30.  
  31.      BBBB REAL(N)
  32.         the right hand side vector.
  33.  
  34.      JJJJOOOOBBBB INTEGER
  35.         = 0         to solve  A*X = B ,
  36.         = nonzero   to solve  TRANS(A)*X = B  where
  37.         TRANS(A)  is the transpose.  On Return
  38.  
  39.      BBBB the solution vector  X . Error Condition
  40.  
  41.      AAAA division by zero will occur if the input factor contains a zero on the
  42.      diagonal.  Technically, this indicates singularity, but it is often
  43.      caused by improper arguments or improper setting of LDA .  It will not
  44.      occur if the subroutines are called correctly and if SGECO has set RCOND
  45.      .GT. 0.0 or SGEFA has set INFO .EQ. 0 .  To compute  INVERSE(A) * C
  46.      where  C  is a matrix with  P  columns
  47.         CALL SGECO(A,LDA,N,IPVT,RCOND,Z)
  48.         IF (RCOND is too small) GO TO ...
  49.         DO 10 J = 1, P
  50.         CALL SGESL(A,LDA,N,IPVT,C(1,J),0) 10 CONTINUE LINPACK.  This version
  51.      dated 08/14/78 .  Cleve Moler, University of New Mexico, Argonne National
  52.      Lab.  Subroutines and Functions BLAS SAXPY,SDOT
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.